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ABSTRACT 

A  technique  that  consists  of  correlation  followed  by  thresholding  is 
widely  used  to  mitigate  multipath  effects,  especially  in  the 
presence  of  non-line-of-sight  (NLoS)  effects,  for  TOA  estimation. 
In  practice,  an  accurate  threshold  that  the  technique  relies  on  is 
difficult  to  obtain,  since  optimizing  the  threshold  requires  prior 
knowledge  of  channel  statistics  and  signal/noise  power,  which  are 
not  always  known  a  priori  nor  are  easy  to  estimate.  In  this  paper, 
we  propose  a  new  thresholding  method  that  does  not  rely  on  such 
prior  knowledge  and  is  robust  against  an  inaccurate  threshold. 

Index  Terms —  non-line-of-sight,  threshold,  TOA,  multipath 

1.  INTRODUCTION 

TOA  or  TDOA  methods  are  important  means  for  location 
estimation,  but  are  often  plagued  by  multipath  propagation, 
especially  in  non-line-of-sight  (NLoS)  multipath  channels  where  it 
is  difficult  to  identify  the  obstructed  LoS  path  for  TOA  estimation, 
since  the  LoS  path  is  usually  not  the  strongest  path.  In  addition, 
noise  causes  spurious  peaks  that  appear  to  be  false  paths  that 
obscure  a  weak  LoS  path.  Based  on  the  fact  that  these  false  paths 
tend  to  he  weaker  than  the  LoS  path,  existing  thresholding 
techniques  remove  the  false  paths  using  a  threshold,  which  is 
optimized  either  on-line  or  off-line.  The  off-line  method  [1]  relies 
on  prior  knowledge  of  the  channel,  such  as  the  distribution  of  the 
LoS  path  gain.  The  on-line  method  [2]  needs  to  estimate  the 
signal-to-noise  ratio  (SNR),  which  is  not  a  trivial  task.  Moreover, 
thresholding  methods  are  not  robust  against  an  inaccurate 
threshold,  and  if  the  threshold  is  inaccurate,  the  possibility  of 
missed  and  false  detection  of  the  LoS  path  can  be  large  (cf.  Fig.  5). 

The  central  idea  of  our  proposed  method  is  that  it  is  unlikely 
that  two  false  paths,  which  are  caused  by  a  high  level  of  noise,  in 
two  independent  TOA  estimates  have  the  same  time  delay;  while 
the  time  delays  of  the  LoS  path  are  likely  to  remain  the  same  across 
different  TOA  estimates.  In  other  words,  the  time  delays  of  the 
false  paths  and  those  of  the  LoS  path  exhibit  different  kinds  of 
statistics.  Based  on  this  property,  we  introduce  an  integer 
threshold,  by  optimizing  which  the  LoS  path  can  he  distinguished 
from  false  paths.  Our  method  is  not  dependent  on  an  accurate 
threshold  and  precise  knowledge  of  the  channel,  and  is  therefore 
robust  and  can  be  applied  to  a  wide  range  of  signals. 
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2.  SIGNAL  MODEL 

We  present  a  signal  model  for  the  reception  of  training 
sequences  (TRS)  in  a  multipath  environment.  The  baseband 
transmitted  signal  s,(t)  for  the  i  th  TRS  is  formed  by  modulating 
symbols  with  the  symbol  waveform  g(t) , 

N+N<- 1 

sid)  =  E  \n-N,)9d  -  nT )  (1) 

n=Ni 

where  the  TRS  of  length  L  is  assumed  to  he  a  pseudo¬ 

random  noise  (PN)  code,  and  T  is  the  symbol  period.  The  i  th 
TRS  starts  from  the  Ni  th  symbol  in  the  data  stream,  where  N,  is 
unknown.  The  data  stream  is  formatted  into  segments  of  lengths 
(Ai+J  —  Nt)T  ,  with  the  same  length  LLTRS  at  the  same  position 
in  each  segment  (Fig.  1).  We  assume  that  the  time  duration  LT  of 
one  TRS  is  much  smaller  than  the  channel  coherence  time,  and 
thus  the  time-varying  fading  channel  gains  remain  unchanged 
within  LT  hut  may  vary  over  different  TRSs.  The  channel 
impulse  response  (CIR)  hit)  corresponding  to  the  z  th  TRS  is 
given  by 

K,-1 

W )  =  %sd  -  To,i)  +  E  ckAf  ~  tm)  (2) 

k= 1 


where  the  fcth  path  is  parameterized  by  a  path  delay  r,  .  and  a 
complex  path  gain  ck  t ,  and  Kp  is  the  number  of  multipaths.  The 
parameters  r0 ,■  and  c0  ■  correspond  to  the  propagation  delay  and 
the  gain  of  the  LoS  path,  respectively.  The  segment  of  the  received 
baseband  signal  corresponding  to  the  i  th  TRS  is  given  by  [3] 


n{t)  =  s,(t)  ®  \{t)  +  zfi) 


Kp- 1 

E  c*,a(* 


k= 0 


+  z,(t) 


(3) 


where  the  channel  noise  zt(t)  is  zero-mean  additive  white 
complex  Gaussian  with  a  constant  power  spectral  density  N0  .  The 
carrier  frequency  is  assumed  to  be  estimated  and  compensated,  so 
that  there  is  no  phase  rotation  term  [3], 

Y  (F+i  -  N,)T 
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Fig.  1  Data  symbol  stream  with  TRSs 
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3.  ALIGNMENT  OF  CORRELATOR  OUTPUTS 

We  cross-correlate  the  received  signal  with  a  locally 
generated  TRS  PN  waveform  via  a  channel  correlator  [3],  The 
locally  generated  PN  waveform  x(t)  is  formed  by  multiplying  the 
PN  code  with  g{t) ,  i.e.,  x(t)  =  En= ob^  —  nT)  .  Assuming 
that  the  PN  code  has  an  ideal  autocorrelation  function,  and  that  the 
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cross-correlation  between  the  PN  code  and  all  other  data  symbols 
is  zero,  the  cross-correlation  can  be  written  as  a  convolution: 

kp-  t 

Vi(t)  =  n(t)  ®  x*{-t)  =  A  J2  ckjp{t  -  N{T  -  Tkii)  +  nit)  (4) 

k= 0 

where  A  =  L  ■  E^bn  |2}  is  a  constant,  the  correlation  function  is 
p(t)  =  g(t)  ®  ( (■)*  denotes  complex  conjugate),  and  the 

Gaussian  stationary  noise  nt (t)  is  colored  with  the  known 
autocorrelation  rn  (r)  =  AN0p(r) . 

We  truncate  and  sample  yit)  to  obtain  a  finite  length  signal 
of  duration  2 (  L  „  is  the  maximum  delay  spread  and  is 
known  a  priori),  which  contains  two  Lmax  portions  before  and 
after  the  strongest  path,  since  the  weak  LoS  path  may  he  located 
before  the  strongest  path  due  to  NLoS  effect.  In  this  way,  the 
truncated  correlator  outputs  {yit)}  are  guaranteed  to  contain  the 
LoS  path,  but  they  may  also  include  false  paths  caused  by  channel 
noise.  The  finite  length  {yAt)}  for  different  time  segments  i  need 
to  he  aligned  for  the  rest  of  the  method  to  work  properly.  There  are 
three  cases  to  consider. 

Case  I:  transmitter  and  receiver  are  stationary.  It  is  seen  from 
(2)  and  (4)  that  the  TOA  of  the  i  th  TRS  is  (AFT  +  r0i)  . 
Although  NT  is  unknown,  in  many  cases  the  TRS  repetition 
period  (AF+1  —  AF)T  is  known.  And  the  LoS  path  propagation 
delay  rQi  is  unchanged  over  different  TRSs,  because  transmitter 
and  receiver  are  stationary.  Therefore  the  (i  +  1)  th  correlator 
output  can  be  aligned  with  that  of  the  i  th  with  a  time  shift  of 

(Ni+i  -Nt)T. 

Case  II:  transmitter  or  receiver  is  moving.  There  is  still  a 
mismatch  between  consecutive  correlator  outputs  after  the  time- 
shifting  in  Case  I,  since  the  LoS  path  propagation  delay  roi 
changes  as  transmitter  or  receiver  moves.  The  mismatch  linearly 
increases  or  decreases  if  the  motion  is  in  constant  speed.  We  can 
add  a  parameter  drQi  /  dt  that  characterizes  the  linearity,  and 
jointly  estimate  r0  i  and  drQ  i  /  dt ,  which  can  he  tracked  via  a 
Kalman  filter.  In  each  of  the  iterations,  the  alignment  is  done  via 
dr0  i  /  dt  that  is  estimated  in  the  previous  iteration.  Due  to  space 
limit,  we  will  not  discuss  this  recursive  approach  further. 

Case  III:  (Ni+l  —  Nf)T  is  unknown.  The  methods  in  Case  I 
and  II  can  not  be  applied  to  this  case.  But  if  we  assume  there  is  a 
strong  correlation  between  consecutive  correlator  outputs,  i.e.,  two 
consecutive  TRSs  are  within  the  channel  coherence  time,  which 
means  (Ni+1  —  Nt)T  is  relatively  small,  a  sequential  alignment 
algorithm  [4]  can  he  applied. 

In  the  remainder  of  the  paper,  we  assume  that  the  alignment 
of  y, (t)  has  been  performed,  which  will  he  used  in  Section  5. 

4.  THEORETICAL  OPTIMAL  THRESHOLD 

In  NLoS  channels,  the  TOA  estimation  approach  in  [1,  2]  is 
to  set  a  threshold  tpth  to  remove  all  the  weak  paths  below  the 
threshold,  and  search  over  the  finite  duration  2 £  of  the 
estimated  CIR  envelope  \h-(t  —  NtT) |  estimated  from  yt(t)  for  the 
first  peak  position,  which  is  supposed  to  be  the  estimated  TOA  of 
the  i  th  TRS. 

However,  if  we  set  the  threshold  ip.h  too  high,  there  is  a 
possibility  of  the  LoS  path  being  removed.  This  probability  of  a 
missed  detection  of  the  LoS  path  can  be  written  as  [2] 

PM  =P*{HNiT  +  T^\<^h}  (5) 


where  |  ^(AFT  +  roi)  is  the  peak  strength  of  the  correlator 
output  corresponding  to  the  LoS  path.  On  the  other  hand,  if  the 
threshold  rith  is  too  low,  we  may  fail  to  remove  a  false  path 
positioned  before  the  LoS  path.  This  will  result  in  a  false  detection. 
As  shown  in  Fig.  4  (h)  (c),  a  false  path  is  a  false  peak  caused  by  a 
high  level  of  noise,  and  the  probability  of  a  false  detection  can  be 
approximated  by  [  1] 


(6) 


where  F  (of  length  Ln  )  indicates  the  noise-only  portion  (or 
“noise  portion”  for  simplicity,  c.f.,  Fig.  4)  of  y^t) .  Eq.  (6)  is  the 
probability  of  the  envelop  of  the  stationary  process  n;(f)  reaching 
above  the  level  in  a  time  duration  L  .  With  a  probability  of 
PF  +  PM  ,  &  false  detection  or  a  missed  detection  causes  large 
errors  (outliers)  of  TOA  estimates.  PF  +  PM  can  be  viewed  as  a 
Bayes  risk,  by  minimizing  which  we  can  determine  the  optimal 
threshold  tl>th  to  minimize  the  number  of  outliers. 

Calculation  of  PM  is  not  an  easy  task,  since  the  distribution 
of  yiNjT  +  roi)|  relies  on  both  signal/noise  power  and  channel 
statistics  (such  as  the  distribution  of  the  LoS  path  gain)  [1,  2],  In 
contrast,  it  is  relatively  easy  to  compute  PF  .  Note  that  since  PF  is 
a  function  of  a  small  time  interval  T ,  it  does  not  simply  depend  on 
the  PDF  of  nt(t)  .  Rather,  it  depends  heavily  on  the  autocorrelation 
of  nit)  which  is  known  (the  channel  noise  power  can  be 
estimated  cf.  the  next  paragraph).  As  shown  in  Fig.  2,  we  compute 
the  PF  by  modeling  the  noise  process  n;(f)  as  a  two-state 
Markov  process,  which  falls  in  state  1  when  jn^f)  >  ipth  for 
some  t\  and  otherwise  is  in  state  0.  This  approach  is  analogous  to 
the  finite-state  Markov  channel  model  [5],  where  the  continuous 
time-varying  channel  is  quantized  into  a  set  of  intervals  (states). 
The  difference  is  that  the  channel  autocorrelation  is  usually 
characterized  by  the  Jakes  filter. 


Fig.  2  Two-state  Markov  process  for  channel  noise 


Next  we  will  determine  the  transition  rates  A  and  p,  of  the 
two-state  Markov  process.  We  define  n^t)  =  ^(f)  +  jr]2(t)  , 
where  Gaussian  processes  ^(t)  and  rj2(t)  are  uncorrelated  with 
identical  autocorrelation  functions,  i.e.,  r  (r)  =  rl;  (r) 
=  (1  /  2 )rn  (t)  .  The  level-crossing  rate  of  | nt(t)\  crossing  the 
level  tj)th  is  given  by  [6] 


LCR(ipth) 


2<Zn 


(7) 


where  (Jq  =  (0)  denotes  the  power  of  ??1(£)  and  r]2(t)  ,  and 
/3  =  —r  (0)  (  irv  ( t )  =  d2rv  (r)  /  dr2  ).  Note  that  a2  can  be 
estimated  by  finding  a  level  p  that  maximizes  the  number  of  level 
crossings  of  y.(t)  in  the  first  half  of  2 L  (consisting  of  mostly 
noise),  since  argma xpLCR(p)  =  o0 .  The  average  length  of  the 
time  intervals  in  which  the  envelop  of  n;(i)  is  below  a  given  level 
tpth  ,  i.e.  the  average  period  of  state  0,  1  /  A  ,  is  given  by  [6] 


1  Pr {\nit)\  <  iPth} 

few  <Jq 

exp 

Ah 

-  i 

A  LCR^J 

V  P 
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Similarly,  the  average  period  of  state  1  is  given  by 

1  /  M  =  pr{jkW|  >  V’tfJ  /  LCR^) . 

The  noise  portion  length  L  is  assumed  to  be  a  fixed  value 
since  aligning  y^t)  makes  Ln  of  different  yt(t)  equal.  Noticing  (6) 
is  the  same  as  that  of  the  two-state  Markov  process  visiting  state  1 
at  least  once  in  the  time  duration  L  ,  Pp  that  is  needed  in  Fig.  5 
is  derived  as  (using  the  properties  of  the  Markov  process  [5]) 


P,  =  1  - 


A  +  p 


(9) 


5.  COARSE  DETECTION  OF  LOS 

As  shown  in  Fig.  5,  PF  +  PM  vs.  the  threshold  i/)th  is  a  U 
shaped  curve,  the  bottom  of  which  corresponds  to  the  minimum 
F F  +  PM  .  Decreasing  L  ,  which  is  equivalent  to  making  a  coarse 
estimation  of  LOS  TOA,  not  only  reduces  the  minimum  PF  +  PM 
but  also  makes  a  wider  range  for  an  acceptable  threshold. 
Intuitively,  directly  summing  the  envelopes  of  aligned  correlator 
outputs,  i.e.  {|  y.jit)  }  ,  will  increase  SNR  and  therefore  will  have  a 
similar  effect.  However,  the  peaks  in  the  noise  portion  of  the  sum 
of  |}  can  be  large  due  to  increased  variance  by  summation. 

In  our  method  we  will  reduce  this  variance  by  quantizing 
\yi(t)\  to  the  two  values  of  the  Markov  process  and  perform  a 
coarse  detection  of  the  LOS  path,  under  the  condition  that  only  the 
noise  power  is  estimated  and  known.  Because  of  this  condition, 
Pf  +  Pm  's  unknown  and  thus  the  optimal  threshold  cannot  be 
computed.  However,  the  transition  rates  A  and  p  of  the  two-state 
Markov  process  are  controllable  by  .  We  choose  a  threshold 
ipfo  that  satisfies  i pth  >  y]'2\g(2)a0  =  1.18<r0  ,  which  makes 
X  /  p  <1  thereby  Pp  small  (cf.  the  factors  of  A  /  p  in  (12)). 
Threshold  also  needs  to  be  small  enough,  so  that  the  LoS  path 
is  unlikely  to  miss,  i.e.,  PM  is  small.  In  our  method  some  rough 
knowledge,  such  as  in  the  worst  case  scenario  PM  is  10%  when 
ipth  is  less  than  4 ag  (computed  from  (5)),  i.e.,  PM  <  10%  for 
ipth  <  4 cr0  ,  is  sufficient  for  the  rest  of  the  method  to  work.  This 
kind  of  rough  knowledge  of  a  tolerable  PM  is  usually  available  in 
practice.  The  choice  of  this  threshold  is  rather  subjective  and  a 
wide  range  of  the  threshold  can  be  accepted. 

Using  this  rough  threshold  i/jth  ,  each  correlator  output  y.(t) , 
which  contains  both  the  noise  portion  and  the  signal  portion  (the 
region  containing  signals,  the  leading  edge  of  which  corresponds 
to  the  LoS  path.  cf.  Fig.  4(d)),  is  quantized  into  2  states  to  generate 
a  two-state  on-off  process  vt(t)  ,  i.e.  vf(t)  =  0  when  the  envelop 
of  i /.(t)  is  below  the  threshold,  i.e.  |j/;(t)  <  i/)th  for  some  t,  and 
otherwise  ir(f)  =  1  for  some  other  t.  The  quantization  of  y.[t) 
into  vit)  effectively  reduces  the  variance  of  the  noise  portion, 
therefore  v,(t)  is  more  robust  than  y.(t)  even  for  a  very  rough 
estimate  of  PM .  Thus,  rather  than  y.(t),  we  sum  vfi)  to  obtain  a 
discrete  valued  process  v(t) : 

M 

v(t)  =  J2vi(t)  o°) 

«=i 

where  (ir(t)}  are  aligned  and  thus  the  noise  portion  of  v(t) 
contains  only  noise.  In  the  viewpoint  of  packet  multiplexing  [5], 
the  noise  portion  of  v(t)  is  the  sum  of  M  identical  two-state  on- 
off  mini-sources,  and  thus  it  can  be  modeled  by  an  (M  +  1)  -state 
Markov  chain  with  state  dependent  transition  rates  (cf.  Fig.  3). 

The  noise  portion  and  the  signal  portion  of  v(t)  show 
different  statistical  characteristics  (cf.  Fig.  4(d)).  In  the  noise 


portion,  v(t)  tends  to  stay  in  states  of  small  numbers,  since  the 
false  paths  (corresponding  to  state  1  of  v.(t))  over  different  v >.(t) 
are  unlikely  to  locate  at  the  same  position  and  therefore  their  sum 
is  small.  By  contrast,  in  the  signal  portion,  with  a  high  probability, 
v(t)  stays  in  states  of  large  numbers,  e.g.  state  M.  Therefore,  we 
define  an  integer  threshold  m.h  ( 0  <  mth  <  M  —  1 )  for  v(t)  ,by 
minimizing  a  new  Bayes  risk  defined  below,  and  the  first  transition 
time  from  state  mth  to  ( mth  +1)  is  a  coarse  detection  of  the  LOS 
path  arrival  time. 

Similar  to  the  analysis  in  Sec.  4,  we  define  new  probabilities 
of  false  and  missed  detection  for  the  discrete  v(t) .  The  probability 
of  detecting  the  LOS  in  each  w  (t) ,  corresponding  to  state  1,  is 
(1  -  pm)  ,  and  fails  with  probability  pm  ■  A  missed  detection  of 
v(t)  occurs  when  the  number  of  events  of  detecting  the  LoS  path 
for  each  v.(t)  is  less  than  ( m.h  +  1)  .  Therefore  the  number  of 
such  events  (the  sum  of  state  1  of  v.(t))  is  binomial  distributed, 
and  the  probability  of  missed  detection  by  v(t)  is  given  by 


P'm=  E 

k= 0 


(1  Pm)*  Pm 


(11) 


A  false  detection  happens,  when  the  first  transition  from 
state  mth  to  (m.,  +  1)  occurs  in  the  noise  portion.  This  is  a  first- 
passage-time  problem  for  finite-state  Markov  chains.  The  closed- 
form  expression  for  this  kind  of  problem  is  not  always  available 
and  further  discussion  is  beyond  the  scope  of  this  paper.  Instead, 
for  the  case  that  M  is  moderately  large,  we  approximate  the  false 
detection  probability  as  (this  derivation  is  omitted  here  due  to 
space  limitation.) 


p' 

rF 


M 


E 

A=ma+i 


'm' 

a]a 

Ml 

k 

(12) 


Similar  to  the  analysis  in  Sec.  4,  we  can  minimize  Pp  +  P'u 
to  obtain  an  optimal  integer  threshold  mth  .  As  shown  in  Fig.  6  , 
the  minimal  P'  +  P[,  is  much  smaller  than  P„  +  P,,  and  the 
threshold  mlh  can  be  easily  chosen  from  a  small  set  of  integers. 


Fig. 3  (M+l)-state  Markov  process  for  noise  portion  of  v(t) 


Our  coarse  detection  method  is  theoretically  more  tractable 
than  directly  summing  {|j/;(f)|}  which  will  involve  calculating  a 
continuous  threshold  that  is  not  trivial  without  knowledge  of  signal 
power  and  channel.  In  contrast,  determining  an  integer  threshold  is 
much  easier  as  shown  above.  The  coarse  detection  is  more  robust 
and  can  detect  the  LoS  TOA  to  within  a  few  symbols  in  simulation. 

6.  REFINEMENT  OF  LOS  TOA 

If  more  precise  TOA  of  the  LoS  path  is  required,  we  can 
refine  the  above  result  by  using  again  the  correlator  output  y.(t) , 
but  only  in  the  vicinity  of  the  coarsely  detected  LoS  time.  Since  we 
are  now  so  close  to  the  actual  LoS  TOA,  there  is  no  longer  any 
issue  of  false  or  missed  detection  -  the  only  issue  now  is  accuracy. 
So  thresholding  is  no  longer  needed  and  therefore  y.{t)  can  now 
be  used.  But  y  (t)  is  also  subject  to  inter-multipath  interferences 
caused  by  the  tails  of  the  correlation  function  p(t)  (cf.(4)).  If  the 
pulse  p(t)  is  relatively  narrow  compared  to  path  delay  differences, 
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e.g.,  the  waveform  g(t)  is  wideband  and  the  paths  are  not  very 
dense,  the  multipath  interferences  between  paths  are  small.  In  such 
a  case,  the  estimated  CIR  h^t  —  NT)  can  be  approximated  by  the 
peaks  of  y^t)  /  (A|p(0)|) ,  and  the  closest  peak  of  ht(t  —  NtT)  to 
the  coarsely  detected  LoS  TOA  corresponds  to  the  refined  TOA. 
Otherwise,  we  need  to  further  process  yXt)  to  obtain  accurate  CIR 
via  e.g.  the  time-domain  approach  [1]  or  the  frequency  domain 
approach  [7],  and  obtain  a  refined  LoS  TOA  by  the  closest  peak  to 
the  coarsely  estimated  LoS  TOA. 

7.  SIMULATION  RESULTS 

The  simulation  channel  is  modeled  as  a  12-path  time- varying 
Rayleigh  fading  channel  with  NLoS  propagations,  and  the  channel 
varies  over  different  TRSs.  Due  to  the  NLoS  effect,  the  LoS  path  is 
4  dB  below  the  strongest  path,  as  shown  in  Fig.  4(a).  Each  time- 
varying  path  gain  is  generated  according  to  the  Jakes  model  with 
the  Doppler  frequency  100  Hz.  The  signal  is  the  ATSC  digital  TV 
signal,  in  which  a  training  sequence  (511  symbol  PN  code) 
periodically  appears  in  data  stream  for  every  24.2  ms.  We  set  the 
oversampling  rate  to  2  /  T  ,  and  the  maximum  delay  spread 
Amx  =  169T  (T  =  0.093 /is  ).  The  CIRs  are  estimated  based  on 
y.(t)  using  the  method  in  [7].  The  SNR  is  defined  as  the  ratio  of 
the  signal  power  of  the  LoS  path  to  the  power  of  the  white 
Gaussian  channel  noise  within  the  signal  bandwidth  (6MHz).  With 
the  known  channel  statistics,  PF  +  PM  is  numerically  evaluated, 
as  shown  in  Fig.  5.  Fig.  6  is  calculated  from  (1 1)  and  (12).  . 

Tab.  1  shows  the  TOA  RMSE  (root  mean  square  error) 
results  of  our  method  (both  coarse  and  refined  TOA)  and  those  of 
the  thresholding  method  [1]  with  an  optimal  threshold  that  is  pre¬ 
calculated  based  on  the  known  SNR  and  the  known  channel 
statistics  (which  are  impractical).  It  is  seen  that  our  TOA  results, 
with  much  less  information,  are  close  to  those  of  [1]  although  the 
threshold  ipth  =  3u0  is  not  optimal.  It  is  also  seen  at  low  SNR  ( 
-KMB),  the  refined  TOA  is  worse  than  the  coarse  TOA.  This  is 
because  there  are  many  false  peaks  caused  by  the  high  level  of 
noise  in  the  vicinity  of  the  LoS,  the  refinement  (Sec.  6)  that  picks 
the  closest  peak  to  the  coarse  LoS  TOA  tends  to  choose  false  peaks 
instead  of  the  peak  corresponding  to  the  actual  LoS  path. 


100  200  300  400  500  600  700  800  900 

Fig.  4  (a)  Delays  and  relative  powers  of  a  12-path  channel.  (b)(c) 
Aligned  correlator  outputs  (d)  Discrete  valued  process  v(t)  for 
M  =  15  ,  ijifh  =  3fj()  and  SNR=  0  dB.  The  integer  threshold 
mth  =  5  ,  which  is  denoted  by  a  horizontal  strait  line,  (x-axes  for 
(a)  (b)  (c)  (d)  are  in  sampling  periods,  i.e.,  T/2.) 


Ln  =  {2,  32,  62,  92,  122,  152}  x  T  and  SNR  =  {-5,  0}  dB 


Fig.  6  Pp  +  P'M  vs.  integer  threshold  mth  for  M  ={15,  21}  (M  is 
the  number  of  TRSs),  PM  =  {5,  6,  7,  8,  9,  10}%  and  =  3 <j(j 
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Tab.  1  TOA  RMSEs  (in  symbol  periods)  vs.  SNR  for  M  =21  (Mis 


the  number  of  TRSs)  and  iplh  =  3<r0  (600  trials). 
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